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Several mechanisms based on neutrino oscillations can explain the observed motions of pulsars if 
the magnetic field in their interiors is of order 10 14 — 10 lj G. 
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y—t ' The proper motions of pulsars ill] present an intriguing astrophysical puzzle. The measured velocities of pulsars 
exceed those of the ordinary stars in the galaxy by at least an order of magnitude. The data suggest that neutron 
stars receive a powerful "kick" at birth. Whatever the cause of the kick, the same mechanism may also explain the 
| rotations of pulsars under some conditions S. 

The origin of the birth velocities is unclear. Born in a supernova explosion, a pulsar may receive a substantial kick 
due to the asymmetries in the collapse, explosion, and the neutrino emission affected by convection Evolution of 
close binary systems may also produce rapidly moving pulsars Q. It was also suggested || that the pulsar may be 
accelerated during the first few months after the supernova explosion by its electromagnetic radiation, the asymmetry 
resulting from the magnetic dipole moment being inclined to the rotation axis and offset from the center of the star. 
t-H ■ Most of these mechanisms, however, have difficulties explaining the magnitudes of pulsar spatial velocities in excess 
of 100 km/s. Although the average pulsar velocity is only a factor of a few higher, there is a substantial population 
OS ' °f pulsars which move faster than 700 km/s, some as fast as 1000 km/s Q|. 

Neutrinos carry away most of the energy, ~ 10 53 erg, of the supernova explosion. A 1% asymmetry in the distribution 
■ of the neutrino momenta is sufficient to explain the pulsar "kicks" . A strong magnetic field inside the neutron star 
could set the preferred direction. However, the neutrino interactions with the magnetic field are hopelessly weak. 

Ordinary electroweak processes || cannot account for the necessary anisotropy of the neutrino emission JtJ . The 
possibility of a cumulative build-up of the asymmetry due to some parity-violating scattering effects has also been 
considered pi. However, in statistical equilibrium, the asymmetry does not build up even if the scattering amplitudes 
^ . are asymmetric MM- Although some net asymmetry develops because of the departure from equilibrium, it is too 
small to explain the pulsar velocities for realistic values of the magnetic field inside the neutron star Jj],[l0| . 

There is a class of mechanisms, however, that can explain the birth velocities of pulsars as long as the magnetic field 
inside a neutron star is 10 14 — 10 15 G. These mechanisms |n|-|l4[ have some common features. First, the conversions 
of some neutrino v into a different type of neutrino, v' , occurs when one of these neutrinos is free-streaming while the 
other one is not. The free-streaming component is out of equilibrium with the rest of the star, which prevents the 
wash-out of the asymmetry. Second, the position of the transition point it affected by the magnetic field. I will review 
two possible explanations, which do not require any exotic neutrino interactions and rely only on the established 
neutrino properties, namely matter-enhanced neutrino oscillations. The additional assumptions about the existence 
of sterile neutrinos |l5[] and the neutrino masses appear plausible from the point of view of particle physics. 

II. PULSAR KICKS FROM NEUTRINO OSCILLATIONS 

As neutrinos pass through matter, they experience an effective potential 

V{v.) = (1) 
V(v e ) = -V(v e ) = V {3Y e -l + 4Y„ e ) (2) 
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V(^r) = -V{^, T ) = V (Y e - 1 + 2Y Ue ) + \-£ j _ (3) 

where Y e {Y Vc ) is the ratio of the number density of electrons (neutrinos) to that of neutrons, B is the magnetic 
field, k is the neutrino momentum, Vq = 10 eV (/3/10 14 gcm~ 3 ). The magnetic field dependent term in equation (j^) 
arises from a one-loop finite-density contribution JlT| , p~8[ | to the self-energy of a neutrino propagating in a magnetized 
medium. An excellent review of the neutrino "refraction" in magnetized medium is found in Ref. Jl8| . 
The condition for resonant pj| ] oscillation i>{ <-* Uj is 

i2 17)^ 



cos 2% + V{vi) = -JL cos 2% + V{ Vj ) (4) 



where fij can be either a neutrino or an anti-neutrino. 

The neutron star can receive a kick if the following two conditions |ll|-|l3|| are satisfied: (1) the adiabatkf] oscillation 
vi <-* Vj occurs at a point inside the i-neutrinosphere but outside the j-neutrinosphere; and (2) the difference [V(^i) — 
V(i^-)] contains a piece that depends on the relative orientation of the magnetic field B and the momentum of the 
outgoing neutrinos, k. If the first condition is satisfied, the effective neutrinosphere of i/j coincides with the surface 
formed by the points of resonance. The second condition ensures that this surface (a "resonance-sphere") is deformed 
by the magnetic field in such a way that it will be further from the center of the star when (k ■ B) > 0, and nearer 
when (k ■ B) < 0. The average momentum carried away by the neutrinos depends on the temperature of the region 
from which they exit. The deeper inside the star, the higher is the temperature during the neutrino cooling phase. 
Therefore, neutrinos coming out in different directions carry momenta which depend on the relative orientation of k 
and B. This causes the asymmetry in the momentum distribution. An 1% asymmetry is sufficient to generate birth 
velocities of pulsars consistent with observation. 

Let us use two different models for the neutrino emission to calculate the kick from the active-sterile and the active 
neutrinos, respectively. As shown in Ref. these two models are in good agreement. 

III. OSCILLATIONS INTO STERILE NEUTRINOS 

Since the sterile neutrinos have a zero-radius neutrinosphere, v s <-> j> MiT oscillations can be the cause of the pulsar 
motions if m(v s ) > m(i>^^ T ). If, on the other hand, m(v s ) < m(i/ M , r ), v s <-> v^ r oscillations can play the same role. 

In the presence of the magnetic field, the condition ([|) is satisfied at different distances r from the center (Fig. |l|), 
depending on the value of the (fc • B) term in (jij). The surface of the resonance is, therefore, 

r(0) = r + 5 cos 4>, (5) 

where cos cp=(k-B)/k and S is determined by the equation (dN n (r)/dr)S » e (3iV e /7r 4 ) 1/3 B. This yields Jl| 



tx z / dr 



where [i e ps (Stt 2 ^) 1 / 3 is the chemical potential of the degenerate (relativistic) electron gas. 

Assuming a black-body radiation luminosity cx T 4 for the effective neutrinosphere, the asymmetry in momentum 
distribution |l^] is 



Afc 4e fn e dT 



k 3-7T 2 V T dN, 



B, (7) 



Non-adiabatic oscillations are discussed in Ref. jit 
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FIG. 1. Sterile neutrinos produced by oscillations are emitted from regions of different temperatures inside a neutron star. 

To calculate the derivative in (Q) , we use the relation between the density and the temperature of a non-relativistic 
Fermi gas. Finally, 



Afc 4e-\/2 
~k ~ 



1/2 
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B = 0.01 
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(8) 



if the neutrino oscillations take place in the core of the neutron star, at density of order 10 14 gcm~ 3 . The neutrino 
oscillations take place at such a high density if one of the neutrinos has mass in the keV range, while the other one 
is much lighter. 

The mixing angle can be very small, because the adiabaticity condition is satisfied if 
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is smaller than the typical scale of the density variations. Thus the oscillations remain adiabatic as long as sin 2 29 > 
10" 8 . 

IV. OSCILLATIONS OF ACTIVE NEUTRINOS 



The active neutrino oscillations can also explain the pulsar kick [pTj . The magnitude of the kick can be calculated 



using a model for neutrino transfer used in the previous section |11|. That is, one can assume that the neutrinos are 
emitted from a "hard" neutrinosphere with temperature T(r) and that their energies are described by the Stefan- 
Boltzmann law. Alternatively, we can use the Eddington model for the atmosphere which was used by Schinder and 
Shapiro to describe the emission of a single neutrino species. One can generalize it to include several types of 
neutrinos J^] 

In the diffusion approximation, the distribution functions / are taken in the form [Ell: 



2 A recent attempt to use the Eddington model for the neutrino transfer failed to produce a correct result because the 
neutrino absorption v e n — > e~p+ was neglected, and also because the different neutrino opacities were assumed to be equal to 
each other. The assumption fcol that the effect of neutrino oscillations can be accounted for in a simplistic model with one 
neutrino species and a deformed core-atmosphere boundary is also incorrect because the temperature profile is determined by 
the emission of six neutrino types, five of which are emitted isotropically. The neutrinos of the sixth flavor, which have an 
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(10) 



where / ei? is the distribution function in equilibrium, Aj denote the respective opacities, m is the column mass density, 
m = J p dx, £ = cosa, and a is the normal angle of the neutrino velocity to the surface. At the surface, one imposes 
the same boundary condition for all the distribution functions, namely 



0, for £ < 0, 
2/ e «, for f > 0. 



(11) 



However, the differences in Aj produce the unequal distributions for different neutrino types. 

Generalizing the discussion of Refs. plfl to include six flavors, three neutrinos and three antincutrinos, one can 
write the energy flux as 



F = 2n E 3 dE / £d£ VY/* + 
Jo J-i i=1 



(12) 



We will assume that Aj = hf\E' 2 /El) 
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FIG. 2. Oscillations i/ T *-» i/ e result in anisotropic emission of r-neutrinos. The upward-going i/ T 's escape from the region 
with a lower temperature than that of the downward-going r-neutrinos. 



We use the expressions for f Vi from equation (J10|). Changing the order of differentiation with respect to m and 
integration over E and £, and using the fact that f eq is isotropic, we arrive at the result similar to that of Ref. j gjj : 
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(13) 



The basic assumption of the model is that flux _F is conserved. In other words, the neutrino absorptions v e n — > e~p + 
are neglected. Since the sum in brackets, as well as the flux F are treated Q] as constants with respect to m, one 
can solve for T 2 : 



anisotropic momentum distribution, cause negligible (down by at least a factor of 6) asymmetry in the temperature profile. 
When the neutrino absorption is included, the Eddington model gives the same result for the kick as the model with "hard 
neutrinospheres" Q , E3 |. 
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Fm +[—F) (14) 



Swapping the two flavors in equation ( |l4[ ) leaves the temperature unchanged in the Eddington approximation. 
Hence, neutrino oscillations do not alter the temperature profile in this approximation. 
We will now include the absorptions of neutrinos. 

Some of the electron neutrinos are absorbed on their passage through the atmosphere thanks to the charged-current 
process 

v e n — > e~p + . (15) 

The cross section for this reaction is a — 1.8 G 2 p E 2 , where E v is the neutrino energy. The total momentum transfered 
to the neutron star by the passing neutrinos depends on the energy. 

Both numerical and analytical calculations show that the muon and tau neutrinos leaving the core have much higher 
mean energies than the electron neutrinos |^2],^3| . Below the point of MSW |l9| resonance the electron neutrinos have 
the mean energies w 10 MeV, while the muon and tau neutrinos have energies « 25 MeV. 

The origin of the kick in this description is that the neutrinos spend more time as energetic electron neutrinos on 
one side of the star than on the other side, hence creating the asymmetry. Although the temperature profile remains 
unchanged in Eddington approximation, the unequal numbers of neutrino absorptions push the star, so that the total 
momentum is conserved. 

Below the resonance E Va < E Vt . Above the resonance, this relation is inverted. The energy deposition into the 
nuclear matter depends on the distance the electron neutrino has traveled with a higher energy. This distance is 
affected by the direction of the magnetic field relative to the neutrino momentum. 

We assume that the resonant conversion v e «-* v T takes place at the point r = ro + 8(cf));5((f>) — <5ocos0. The 
position of the resonance depends on the magnetic field B inside the star |ll[] : 

e f i e B /dN e 

where N e w Y e N n is the electron density and fx e is the electron chemical potential. 

Below the resonance the r neutrinos are more energetic than the electron neutrinos. The oscillations exchange 
the neutrino flavors, so that above the resonance the electron neutrinos are more energetic than the r neutrinos. 
The number of neutrino absorptions in the layer of thickness 2S((f>) around ro depends on the angle <f> between the 
neutrino momentum and the direction of the magnetic field. Each occurrence of the neutrino absorption transfers the 
momentum E Ue to the star. The difference in the numbers of collisions per electron neutrino between the directions 
ef> and n + <j> is 

Ak e /E Ue = 2 S(4>) N n [<t(£i) - a(E 2 )} (17) 

Y, 7T 2 



= 1-8 G\ [El - E 2 2 ] £ — h Nc cos <j>, (18) 



where IiN e = [c?(ln Ne)/dr\~ l . 

We use Y e f» 0.1, E\ « 25 MeV, E 2 » 10 MeV, [i e w 50 MeV, and h,N e « 6 km. After integrating over angles 
and taking into account that only one neutrino species undergoes the conversion, we obtain the final result for the 
asymmetry in the momentum deposited by the neutrinos: 

f = »-»'^ < M > 

which agrees with the estimates^ Qj2^] that use a different model for the neutrino emission. 



3 We note in passing that we estimated the kick in Refs. Jll|jl2| assuming fi e m const. A different approximation, Y e ~ const, 



gives a somewhat higher prediction for the magnitude of the magnetic field |m| 
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Neutrinos also lose energy by scattering off the electrons. Since the electrons are degenerate, the final-state electron 
must have energy greater than \i e . Therefore, electron neutrinos lose from 0.2 to 0.5 of their energy per collision in 
the neutrino-electron scattering. However, since N e <C N n , this process can be neglected. 

One may worry whether the asymmetric absorption can produce some back-reaction and change the temperature 
distribution inside the star altering our result (|l9|). If such effect exists, it is beyond the scope of Eddington approxi- 
mation, as is clear from equation (|l4|). The only effect of the asymmetric absorption is to make the star itself move, 
in accordance with the momentum conservation. This is the origin of the kick @. 

Of course, in reality the back-reaction is not exactly zero. The most serious drawback of Eddington model, pointed 
out in Ref. pl|], is that diffusion approximation breaks down in the region of interest, where the neutrinos are weakly 
interacting. Another problem has to do with inclusion of neutrino absorptions and neutrino oscillations | pl[ . However, 
to the extent we believe this approximation, the pulsar kick is given by equation (|l^) . 



V. CONCLUSION 



The neutrino oscillations can explain the motions of pulsars. Although many alternatives have been proposed, all 
of them fail to explain the large magnitudes of the pulsar velocities. 

If the pulsar kick velocities are due to v e <-> i/„ T conversions, one of the neutrinos must have mass ~ 100 eV 
(assuming small mixing) and must decay on the cosmological time scales not to overdose the Universe This 
has profound implications for particle physics hinting at the existence of Majorons |25[ or other physics beyond the 
Standard Model that can facilitate the neutrino decay. 

If the active-to-sterile neutrino oscillations |l2j are responsible for pulsar velocities, the prediction for the sterile 
neutrino to have a mass of several keV is not in contradiction with any of the present bounds. In fact, the ~keV mass 
sterile neutrino has been proposed as a dark-matter candidate |2(| . 

Some other explanations |l4| ] that utilize new hypothetical neutrino properties, but use a similar mechanism for 
generating the asymmetry, can also explain large pulsar velocities. 



VI. ACKNOWLEDGEMENTS 



The author thanks E. S. Phinney and G. Segre for many interesting and stimulating discussions. This work was 
supported in part by the US Department of Energy grant DE-FG03-91ER40662. 



[1] See, e. g., J. E. Gunn and J. P. Ostriker, Astrophys. J. 160, 979 (1970); A. G. Lyne, M. J. Salter, Mon. Not. R. astr. Soc. 

261, 113 (1993); A. G. Lyne and D. R. Lorimer, Nature 369 (1994) 127; B. M. S. Hansen and E. S. Phinney, Mon. Not. 

R. Astron. Soc. 291, 569 (1997); J. M. Cordes and D. F. Chernoff, Astrophys. J. 505, 315 (1998); C. Fryer, A. Burrows, 

and W. Benz, Astrophys. J. 496, 333 (1998). 
[2] H. Spruit and E. S. Phinney, Nature 393, 139 (1998). 

[3] I. S. Shklovskii, Sov. Astr., 13, 562 (1970) [Astr. Zh. 46, 715 (1970)]; S. E. Woosley, in The origin and Evolution of Neutron 
Stars, proceedings of the 125th Symposium of the International Astronomical Union, D. Reidel Pub. Co., Dordrecht, 1987; 
H.-Th. Janka and E. Miiller, Astron. and Astrophys. 290, 496 (1994); A. Burrows and J. Hayes, Phys. Rev. Lett. 76, 352 
(1996). 

[4] J. R. Gott, J. E. Gunn and J. P. Ostriker, Astrophys. J. Lett. 160, L91 (1970). 
[5] E. R. Harrison and E. P. Tademaru, Astrophys. J. 201, 447 (1975). 

[6] N. N. Chugai, Pis'ma Astron. Zh. 10, 210 (1984); O. F. Dorofeev, V. N. Rodionov and I. M. Ternov, Sov. Astron. Lett. 

11, 123 (1985). 
[7] A. Vilenkin, Astrophys. J. 451, 700 (1995). 



G 



[8] C.J. Horowitz and G. Li, Phys. Rev. Lett. 80, 3694 (1998); D. Lai and Y.-Z. Qian, (a.stro-ph/9802344| Astrophys. J. 495, 



L103 (1998); ibid. 501, L155 (1998); H.-Th. Janka, in Proceedings of "Neutrino Astrophysics", ed. M. Altmann et al., 



[10] 

[ii] 

[12] 
[13] 
[14] 



[15] 
[16] 
[17] 



[18] 
[19] 

[20] 
[21] 
[22] 
[23] 

[24] 
[25] 
[26] 



Tech. Univ. Miinchen, Garching, 1998 ( |atro-ph/9801320| ). 

A. Kusenko, G. Segr e and A. Vilenkin , Phys. Lett. B 43 7, 359 (1998). 

P. Arras and D. Lai, |astro-ph/9806285| ; |astro-ph/9811371 . 

A. Kusenko and G. Segre, Phys. Rev. Lett. 77, 4872 (1996). 

A. Kusenko and G. Segre, Phys. Lett. B 396, 197 (1997). 

A. Kusenko and G. Segre, Phys. Rev. D59, 061302 (1999). 

M. B. Voloshin, Phys. Lett. B 209 (1988) 360; E .Kh. Akhmedov, A. Lanza and D.W. Sciama, Phys. Rev. D 56, 6117 

(1997) ; D. Grasso, FTUV-98-13, |astro-ph/9802060| ; D. Grasso, H. Nunokawa, and J.W.F. Valle, Phys. Rev. Lett. 81, 2412 

(1998) ; R. Horvat, Mod. P hys. Lett. A13, 
See, e.g. R. N. Mohapatra, [hep-ph/9903261; 



_2379 (1998). 

M. Patel, these Proceedings. 
C. W. Kim, J. D. Kim, and J. Song, Phys. Lett. B 419 (1998) 279. 

V. B. Semikoz and J. W. F. Valle, Nucl. Phys. B 425, 651 (1994); ibid. 485, 545 (1997); S. Esposito and G. Capone, 
Z. Phys. C70 (1996) 55; J. F. Nieves and P. B. Pal, Phys. Rev. D40 1693 (1989); J. C. D'Olivo, J. F. Nieves and P. B. Pal, 
ibid., 3679 (1989); J. C. D'Olivo, J. F. Nieves and P. B. Pal, Phys. Rev. Lett., 64, 1088 (1990). 
H. Nunokawa, V. B. Semikoz, A. Yu. Smirnov, and J. W. F. Valle, Nucl. Phys. B 501, 17 (1997). 

S. P. Mikheev and A. Yu. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 913 (1985)]; L. Wolfenstein, Phys. 
Rev. D 17, 2369 (1978). 

H.-Th. Janka and G. G. Raffelt, Phys. Rev. D59 023005 (1999). 
P.J. Schinder and S.L. Shapiro, Astrophys. J. 259, 311 (1982). 

H. Suzuki, in: Physics and Astrophysics of Neutrinos, ed. by M. Fukugita and A. Suzuki (Springer- Verlag, Tokyo, 1994). 
R.W. Mayle, J.R. Wilson and D.N. Schramm, Astrophys. J. 318, 288 (1987); S.W. Bruenn, Phys. Rev. Lett. 59, 938 
(1987). 

Y.-Z. Qian, Phys. Rev. Lett. 79, 2750 (1997 ]; A. Kusenko an d G. Segre, Phys. Rev. Lett. 79 2751 (1997). 
For a recent review, see, e. g., R. D. Peccei, iep-ph/9807516. 

S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72, 17 (1994); G. M. Fuller and X. Shi, |astro-ph/9810076| ; M. Patel, 
these Procedings. 



7 



